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1 Abstract 


An investigation °f flow induced vibration due to coupling between the fluid flow and the propellants lines 
(1A)A and HP-1) was performed. Various flow rate conditions were studied to check whether flow induced 
vibration was possible due to vortex shedding in both valves and pipe lines. Resonance test was conducted 
for all segments of the LOX-feedline for the preburner under test. In addition, critical values of frequency 
an ve ocity are calculated using a mass damping model. A simple chart characterizing the relation between 
frequency and velocity is developed for each component; i.e. propellant lines, valves and flow meters. It was 
found that flow induced vibration occurs for various segments with flow rates of 113 Ib/s, 275 Ib/s and 40 
lb/s. Even more interesting using critical conditions for buckling, it was found that the valve or pipe mav 
collapse for a flow rate of 275 lb/s and valve height of 10% of pipe diameter. Furthermore, two models for 
the acoustic pressure acting on the segments particularly for the valve are proposed. 


2 INTRODUCTION 


Fiow induced vibration occurs when the natural frequency of the line transporting the propellant and the 
uid flow are the same. This matching of the two frequencies yields to a phenomenon known as res onan ce 
This behavior in many cases yield to a collapse of the structure; such as in engine components and propellant 
lines. In order to understand whether re-design or additional testing is required critical values for the fluid 
frequency and flow rate will be calculated for each segment of the E2-LOX line. 

There are various type of phenomena that may induce vibration on components; vortex shedding, turbulence 
water hammer, acoustic among others. Figure 1 shows a flow chart with other cases. However, we will focus 
on vortex shedding studies and the water hammer problem. Vortex shedding occurs when the flow past 
over an obstacle such as cylinder, sphere or any other object is disturbed; resulting in vortices behind* the 
cylinder. These vortices are moving downstream of the flow at a frequency, f s and when this frequency 
is equal to the natural frequency of the structure, vibration is induced by the fluid flow. Water hammer 
occurs normally during the opening or closing of a valve, and it generates an acoustic wave that propagates 
upstream mid downstream of the valve. This transient phenomenon shows up as a loud noise coming out of 
the pipe. This is what you have heard sometimes when you open the water faucet in the morning. Fluid 
ow through valves, bends and orifices generates turbulence as the flow passes through the obstacle. This in 
ino«^ radlates acoustlc waves (of velocity Ua and pressure Pa) upstream and downstream the Valve (Blake 
1996). So as the area of the valves and flow meters changes so the acoustic waves, and when the frequency 
o these waves matches the natural frequency of the pipe, vibration of the components is observed. This is 
because the waves have an acoustic pressure which acts against the surface of the pipe. 

Consequently, the fluid flow and the solid surface are coupled through the forces exerted on the wall by the 
fluid flow. The fluid forces cause the structure to deform, and as the structure deforms it then produces 
c anges in the flow. As a result, a feedback between the structure and flow occurs; action and reaction, 
igure 3 shows a diagram of this action-reaction feedback loop between the fluid flow and structure. Because 
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Figure 1: Flow Induced Vibration Cases. 

of the interaction between the fluid flow and the solid surface the equations of motions describing the 

t yn ?^ CS are j . C ° Up ed ' , This makes the P roblem m °re challenging, and even worse when the flow is turbulent. 
In addition, it means that we must solve the Navier-Stokes equation and the mass damping equation for the 
solid surface simultaneously with their corresponding boundary conditions. 


3 The scope 

As mentioned in previous section flow induced vibration is a result of the coupling of the frequency generated 
by the vortex sheddmg from fluid flow and the natural frequency of the feedlines of the propellants The 
analysis of this investigation is divided into three phases: 

• Phase 1- Acoustic Model: 

This phase consists of modeling the acoustic pressure acting on the structure. The Navier-Stokes 
equation, mass-spring model and structure equation will be used in order to couple the interaction 
e ween t_e -Ui_ „ow and the solid surface. Although the complete solution will not be given at 

this moment, future investigation will be carried out to study the accuracy of the model with actual 
experimental data. 

• Phase 2- Vortex Sheddmg: 

Li this phase we will determine the range of frequencies between the vortices and the natural frequencies 
of the feedlmes. Then we Will determined if flow induced vibration will occur due to coupling between 
the fluid flow vortices and the feedlines. The calculations of this phase will then tell us the critical 
frequency and now rate at which such phenomena occurs. 

• Phase 3- Water hammer 

In this phase flow induced vibration by water hammer phenomena will be studied. Because this 
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Figure 2: FeedbacI loop in Flow Induced Vibration, 


phenomena occurs during the opening or closing of valve, special attention will be placed on valves. 
Determining the critical velocity and frequencies at which resonance occurs is imperative in this phase. 
In addition, buckling of the pipe or collapse of the valve will be carried out. 


4 The Equations of Motion 


It is clear now that the equations of motion for an incompressible flow can be written as, 


and the Navier-Stokes equation, 


dUi 

dx 


= 0 


(i) 


^ + JLr dU - 

dl j dxj pdx'i dxj dxj 


m 


where i and j vary from 1 to 3, thus yield the three components in the x,y and z plane of the velocity field. 
The first term in the left hand side represents the local change of velocity and the second term the transport 
of momentum. The first term from the right hand side is the pressure gradient and the second term is the 
viscous stress term. This is the well known Navier-Stokes equation for incompressible flow. 


The equation describing the deformation of the solid material is obtained from Newton’s second law and is 
given by. 


El^r + pAU 2d2y 


-d 4 y 

dx* 


w +2 ” AU ^ 


+ I™ = ffat) 


QfQ — J *7 (3) 

The first term from the left represents the elastic force experienced by the material as it deforms. Because 
now the material has been deformed, as the fluid flow travels inside the pipe it will experience two new forces 
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arising from the centrifugal acceleration (second term) and coriolis forces (third term). The very last term 
from the left hand side is the acceleration or inertia term, details about this equation is given by Au-Yang 
(2001). Finally but not least, the first term from the right hand side represents the external forcing term 
due to the pressure fluctuations acting on the walls arising from the flow induced vibration. The advantage 
of using this equation is that as the pipe bends or in cases of elbows failure of the material due to buckling 
can be studied. 

It is very clear that these equations are complicated to solve and must be solved simultaneously. Moreover, 
it is important to point out that this equation must be written with their proper bo un dary conditions for 
each element of the pipe including the valves, elbows, and flow meters. However, a simple mass-spring model 
for the system may allow us to provide a simple description of the system. 


4.1 The Mass-Spring Model 

Instead of trying to solve the Navier-Stokes equation and the equation for the solid surface it is possible 
to write down a simple set of equations for the acoustic pressure, Pa, and a mass spring model which 
will essentially provide information about the external forcing term due to flow induced vibration and the 
amplitudes of the wave for various modes. Thus, 

O^v 3v 

(m + nn a )— - + (C s -f- Ci,)— -f- ky — f(t , x) (4) 

The first term is the inertia term or acceleration and it includes the mass added term, m a , which the pipe 
experienced in bends. The second term is similar to the shear stress terms of the Navier-Stokes equation. 
Basically, this term is the damping term due to friction forces at the wall. The last term is the stiffens 
of the material, spring. Finally the only term in the right hand side is the forcing term and comes from 
the external forces due to flow induced vibration. This term needs a close form, and will be related to the 
acoustic velocity, u ay of the wave traveling upstream and downstream of a given point. 

As the fluid flow travels through valves, bends and orifices, it generates turbulence and radiates energy 
upstream and downstream of the area reduced. This acoustic waves reflected by changes in area creates a 
series of standing waves which in turns exert forces on the bends, valves and other surfaces and result in 
vibration. If the acoustics source possesses sufficient energy and the acoustic natural frequencies of the fluid 
flow (i.e. vortex shedding) coincide with the structural frequency then large acoustic motion can occur. In 
fact , I think that this wave propagates upstream and changes the thermodynamics properties of the system, 
resulting in changes of thermodynamic equilibrium points. Figure 3 shows a diagram of the waves generated 
at the valve and propagates upstream toward the tank, causing the conditions at the tank to vary from its 
equilibrium state. 


4.2 The Acoustic Velocity Model 

Based on the above information it makes sense to have a model for the forcing term, /(£,#), which is a 
function of the acoustic wave speed, u a , and the sound speed, c. Thus, 

f(x,t) = P a /A = (Ua j n j )?j (5) 

where rij is the unit vector in the direction of the acoustic velocity, Udj, and the dot product yields to a 
scalar. Notice that an equation for the acoustic wave velocity is required in order to know the forcing term, 
/ (£> or the acoustic pressure. Using the Navier-Stokes equations and the continuity equation it is possible 

to construct an equation for the acoustic wave or for the pressure wave. Assuming that the actual flow has 
the contribution of the mean flow, U Oj, and the acoustic wave, Uaj. This way of decomposing the flow is 
similar to the Reynolds decomposition or linear instability analysis. Thus, 

Uj — Uoj + Uaj ( 6 ) 
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Figure 3: Acoustic Wave in Pipes 


where Uo, » Uaj. Thus, substituting equation 6 into the Navier-Stokes equation and using the continuity 
equation we get, 

dUa dUjUa Id Pa d . dUa, 

dt dxj p dxi + dxj V dxj ^ 

but, the waves are nearly inviscid so the last term can be dropped. Furthermore, this equation can be written 
in terms of the pressure waves and simplified since the mach number. M„ is less than 1 for 
flow. Thus, ' " 

dPa dPa n 

-ZT+C— = 0 fsi 


where c is the sound speed and since we are only interested in one direction (i.e. x\ = x) we then re-write 
the above equation as, 

8Pa dPa „ 

’W +c ~dT = 0 ( 9 ) 

This equation can be solved by using separation of variables and its solution is the particular solution of 
equation 4 and equation 3. It is then necessary to construct a set of equations with its boundary conditions 
for each segment of the feedlines. 


4.3 The Lift and Drag Model 

The goal is to have another expression for the forcing term, f(x,t), that will enable us to study how far the 
acoustic wave will propagate upstream and downstream from a given point. The rationale for the following 
analysis is based on the fact that as the fluid flow passes through a valve or any other obstacle vortex 
structures are form. They in turn produce a lift force and drag force on the object. Thus, the equation for 
the lift force acting on the object is then given by, 

pU 2 

F l = —~DC L sin{2'Kf L t) ( 10 ) 
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where U a is the acoustic velocity to be determined from the Navier-Stokes equations, D is the diameter of 
the object, Cjp is the lift coefficient, and /x, is the frequency at which the lift force acts on the object and 
is equal to the frequency of the vortex shedding, f $ . Moreover, the lift force acts perpendicular to the flow 
direction. Similarly, the object will experience a drag force, given as, 

pU 2 

Fd = ~Y' D( ^Dsin{2'Kf D t) ( 11 ) 

where fjo is the drag frequency of this component given as 2/ 5 , and Cjp is the lift coefficient. Thus, having 
these two forces acting on the object we can then say that total force acting on the object by the acoustic 
wave is, 

/OM) =F l +F d (12) 

now we only need to have a close form for the acoustic velocity. In previous section we relate the acoustic 
velocity with the acoustic pressure (i.e. U a = P a /pc ), and assuming that the actual flow is composed of an 
undisturbed base flow and the wave disturbance due to vortex shedding. Moreover, it is im portant to assume 
in the Navier-Stokes equations that the acoustic wave is inviscid and that the fluid flow is incompressible. 
Therefore, after some simplification and manipulating the equations, the equation for the acoustic velocity 
is given by, 


dUa dUa 


-0 (13) 

where U 0 is the bulk velocity. Using this equation into any of the structure equation (i.e, eqn. 3 or eqn. 4) 
with corresponding boundary conditions will give us the response of the structure due to the acoustic wave 
acting on the surface. 


5 Test Facility and Input uata 


The study of flow induced vibration will be carried out on the E2 facility for the LOX feedline. The schematic 
of this facility is shown in figure 4. Figure 5 shows the diagram for the E2 RF-1 facility. The study will 
evaluate each component shown in the diagrams (i.e. pipes and valves) for various flow rates. Input data 
for the E2 LOX is shown in figure 6 and figure 7. For the LOX studies, flow rates of 40 Ibm/s, 113 ibm/s 
and a maximum of 275 lbm/s were analyzed and for each flow rate studies of resonance for the valve at 
various percentage of valve opening are considered (water-hammer). Some of the empirical calculations and 
approximations used for the model are given in this subsequent sections. 


6 Results and Discussion 


The results of this investigation consists of determining the natural frequency of each segment of the LOX-line 
and the frequency of the vortex sheddings. Studies using the frequency of the simple mass-spring equation will 
be reported in this section. In addition, the frequencies using the material deformation equation, particularly 
buckling of the pipe due to high fluid velocity are computed. The studies considered three cases of flow rates: 
113 Ib/s, 275 Ib/s and 40 lb/ s. 

The frequency from the mass spring model of equation 4 is given by, 


27 rVm 


(14) 


where K represents the stiffness of the pipe and m the mass of the pipe and fluid. The values for K where 
obtained from the finite element analysis using Algor. The values for the frequency of the vortex shedding 
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E2 - LOX LINE 



were determined from the following equation, 



(15) 


where S is the Strouhal number, U is the flow velocity and <U the inner pipe diameter. For the high Reynolds 
number range, 5.43 x 10 5 < Ren < 6.06 x 10 6 , of this research a Strouhal number of about 0.41 was 
appropriate (Blevins, 1994). 


6.1 Vortex Shedding Studies Using Simple Mass-Spring Analysis 

Table 8 shows the flow induced vibration results due to vortex shedding. The first col umn from the left 
refers to the segment number of the LOX line, and the second column is the natural frequency of the pipe 
using equation 14. The third column from the left is the frequency of the vortex shedding (eqn. 15) for a 
mass flow rate of 1131b/s. The fourth column is the frequency for a maximum flow rate of 275 lb/s, and 
the last column is the frequency of the vortex shedding ( eqn. 15) for a flow rate of 40 lb/s. The cells in 
greens indicates that flow induced vibration will occur for a given segment. Clearly most of the flow induced 
vibration due to vortex shedding will occur for a flow rate of 275 lb/s. Notice the vibration for 113 lb/s and 
40 lb/s wifi be expected for segments 3 to 6, while for a flow rate of 275 Ib/s (the maximum for this test, 

RS76) segments 2, 5 and 7 to 9 flow induced vibration will occur. The results from this table are plotted in 
figure 9. 


6.2 The Effect of Valve on Flow Induced Vibration: Wat er- fin m nier 


The focus of this section was to investigate how closing or opening the valve induced vibration on the 
segments of the LOX-lines. Due to the lack of time only the valve located at segment 1 was analyzed. 
Table 10 shows the frequencies where flow induced vibration occurs. The first column is the segment for the 
LOX system. The second one refers to the % of the valve open. The third column is the diam eter of the 


7 




E2 - RP-1 LINE SCHEMATIC 


Key 

^ Rigid Support - All Movement Constrained 


a 

i 


Non-Rigid (Sliding) Support — Axial Movement and 
Movement in Direction of Arrow is Not Constrained 

Line Segment Number 


RS76 CONFIGURATION 



Figure 5: E2-RP1 feedline diagram. 


valve in inches. The fourth column is the natural frequency of the pipe for each segment. The last three 
columns then refer to the vortex shedding frequency of the flow for mass flow rate of 113 lb/s, 275 ib/s and 

40 lb/s consecutively as it passes through the valve. 

Figure 11 shows the results of these same effects. It is clear that the flow rate of 40 Ib/s indeed induced 
£° r the first few segments (i.e. 1 to 6). For higher flow rate such as 113 Ib/s and 2.7.5. .lb/s vibration 
occurs at the segments near the end (i.e. 7 to 9). Notice that the frequencies are not exactly the same as 
the natural frequency of the pipes but close enough that resonance will take place. The fact that resonance 
occurs at 40 Ibm/s may be due to the fact that as the flow is restrained in the valves its local Reynolds 
number increases to a new flow regime, which could be in transition to turbulence. As a result,” other 
phenomena maybe superimpose on acoustic Wave, resulting into a larger disturbance propagating upstream 
and downstream in the LOX-feedlines. 


6.3 Buckling of Pipe due to High Velocity 

In this section we will investigate buckling effects of the pipe as the speed through the valve is increased. 
For this study we will use equation 3 (i.e. the deformation of pipe due to high speeds). Thus, the natural 



Segment 

Description 

O.D. (D 0 )(m.) 

I (in 4 ) 

I.D. (Dt)(m) 

Length (in) 

Pipe Area (in 2 ) 

1 

6” pipe-ACK4 

6.625 

79.155 

4.209 

101.0 

20.558 

2 

6” pipe-ACK4 

6.625 

79.155 

4.209 

177.0 

20.558 

3 

6” pipe-ACK4 

6.625 

79.155 

4.209 

115.5 

20.558 

4 

6" pipe-ACK4 

6.625 

79.155 

4.209 


20.558 

5 

6” p.ipe-ACK4 

6.625 

79.155 

4.209 

29.50 

20.558 

6 

6” pipe-ACK4 

6.625 

79.155 

4.209 

278.3 

20.558 

7 

6” pipe-ACK4 

6.625 

79.155 

4.209 

89.00 


8 

6” pipe-ACK4 

6.625 

79.155 

4.209 

56.50 

20.558 

9a 

6” pipe-ACK4 

6.625 

79.155 

4.209 

10.50 

20.558 

9 

4” pipe-ACK4 

4.500 

17.893 

2.598 

67.40 

■BHI 


Figure 6: E2-L0X Pipe Data: The letter I in the table represents the moment of inertia of the pipe. 


Segment 

Fluid Weight (lb.) 

Pipe Weight (lb.) 


Total Weight (lb) 

Fluid Area (in 2 ) 

1 

57.8 

596 

- 

653.8 

13.914 

2 

101.2 

144.4 

- 

1146 

13.914 

3 

66.0 

681.5 

300 a 

1048 

13.914 

4 

118 

1221.5 

- 

1340 

13.914 

5 

16-9 

r 174.1 


191 

13.914 

6 

1 159 

1642.2 

400 F" 

2201 

13.914 

7 

50.9 

525.2 


576 

13.914 

8 

32.3 

333.4 


366 

13.914 

9a 

6.0 

62 

- 

366 

13.914 

9 

14.7 

205 

800 c 

8288 rf 

5.301 


Figure i . E2-LOX Fluid Data: a-weight of fill line, b-weight of flow meter and filter , c- weight of run valves 
and rigid frame, d-weight includes segment 9a. 


frequency of the structure is given by 


7T I El 
2L 2 y rn p + rrif 


(16) 


where E is the modulus of elasticity I the moment of inertia, m p the mass of the pipe and rrif the mass of 
the fluid. Similarly Au-Yang (2001) determined a critical velocity where buckling of the pipe will take place 
and is given as, 




(17) 


A fluid flow reaching this critical velocity buckling of the pipe should be expected (AU-Yang, 2001). Table 12 
shows the tabulated values for the natural frequency of the pipe using the equation 16. These results represent 
the vortex shedding of the fluid flow for fully open valve case. 


As before, the first column from the left side refers to the segment of the E2-LOX feedline. The second 
column is the natural frequency of the pipe using the equation 16 for material deformation. The next 
three columns are the vortex shedding frequency (eqn 15) for mass flow rate of 113 lb/s, 275 Ib/s, and 40 
Ib/s. The last column from the left is the length of the pipe. The entries in blue mean that flow induced 
vibration occurs for this length. The important point about this table is that only lengths greater than 4 
meters experienced flow induced vibration. Also, notice that essentially flow induced vibration due to vortex 
shedding occurs at flow rate of 275 lb/s which is the maximum expected for this test at the E2 (RS76). 

Table 13 shows the critical velocity of equation 17 for the fully open valve. Notice that no buckling due to 
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high speed will ever occur for any of the flow rates considered. The velocities shown in this table have units 
of m/s. 

However, the results for various valve height are shown in table 14. Notice that when the valve is open 10% 
a velocity of 37m/s is achieved by the flow passing through the valve, and its value exceed the critical speed 
of 21 m/s. Therefore, for a mass flow rate of 2751b/s the segment 1 or the valve will collapse due to high flow 
velocity. The other items mark in purple color does not exceed the critical velocity but, it is close enough 
(Le. 15 m/s for segment 1 at 113 lb/s and 9.47 m/s for segment 2 at 275 lb/s) so that close monitoring of the 
system should be considered. Notice that for segment 2 the critical velocity is 12.07 m/s and the velocity of 
the flow through the valve is 9.47 m/s. For the other segments values as close as these ones are not observed. 

It is very important to understand that the solutions shown here are assuming that a valve of spherical shape 
is used to constrain the flow. Moreover, the nature of the vortex shedding and resonance will depend on the 
geometrical and other variables of the valve. For instance, figure 15 shows some profiles from Castillo at el. 
(2000, 2001) investigating the effect of initial conditions on the downstream flow. Notice that various curves 
for the mean and Reynolds stresses profiles exist depending on the shape, location and wind tunnels speed of 
the tripping device; which can be viewed as different shape and velocities through the valve. The top figure 
is the mean velocity deficit profiles for adverse pressure gradient, favorable and zero pressure gradient flow. 
The interesting observation is that even tough some of the profiles have the same range of local Reynolds 
number and same strength of pressure gradient the profiles are different. This difference observed in all of 
the profiles is due to the fact that different initial conditions were used in each experiment, and as a result 
the profiles of the turbulent flow develop differently depending on those conditions. The profiles at the 
center show the growth of the boundary layer for the adverse pressure gradient profiles of the first figure. 
Again notice the remarkable difference among the profiles. This may explain why the turbulence models 
used under-predict the growth rate of plumes. Finally, the profiles at the bottom represent the Reynolds 
stress profiles in the y-direction for zero pressure gradient data of Castillo and Johansson (2000). Notice 
that this is one of the terms that arise due to turbulence, and a similar effect will be expected in the acoustic 
velocity as the flow is block by the valve or flow meters. 


7 Conclusion 


It was found that flow induced vibration due to vortex shedding do occur for the LOX line of the E2 test 
facility. The analysis for various valve height (water-hammer) and for various flow rate (Le. 113 lb/s, 275 
lb/ s and 40 lb/ s) was performed. Furthermore, for a maximum flow rate of 275 Ib/s bending of the segment 1 
and possibly collapse of the valve may occur. The resonance effect occurring between the pipe and fluid flow 
produced an acoustic wave which propagates upstream and downstream of the test valve. More interesting 
this phenomena may explain the drastic change of the thermodynamic properties of the LOX tank. As 
the flow is restrained by the valve the fluid flow speed increases, and the vortex shedding behind the valve 
increases also. 

Moreover, simple models for the acoustic pressure acting on the pipe surface was proposed. One of the 
approach relates the acoustic pressure with the acoustics velocity. This velocity is obtained from, the Navier- 
S tokos equation. The second approach involved the lift force and drag force experienced by the valve due 
to vortex shedding. They both enter the mass-spring model as a forcing term and its solution describes the 
interaction between the solid surface and fluid flow. 

In summary, further research is required to model the other valve, elbows and flow meter. Testing of the 
models for each segment and component of the LOX line is required. Also, similar analysis is necessary 
for the RPl feedline, particularly of the valves. Verification of this results will be carried usingexperimenta! 
techniques developed by NASA Stennis. Using a software, FREDOM, it will be possible to test the theoretical 
models and results presented in this study. Although some approximations were necessary in order to get 
some understanding of the problem, the present analysis lead to better understanding of the current situation. 
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